Electrical measurement of a two-electron spin state in a double quantum dot 
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We propose a scheme for electrical measurement of two-electron spin states in a semiconductor 
double quantum dot. We calculated the adiabatic charge transfer when surface gates are modulated 
in time. Because of spin-orbit coupling in the semiconductor, spatial displacement of the electrons 
causes a total spin rotation. It follows that the expectation value of the transferred charge reflects 
the relative phase as well as the total spin population of a prepared singlet-triplet superposition 
state. The precise detection of the charge transfer serves to identify the quantum superposition. 
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PACS numbers: 73.21.La 03.67.-a 71.70.Ej 

Electron spins confined in semiconductor quantum 
dots provide intrinsic quantum bits (qubits) In 
2002, an efficient framework for universal quantum com- 
putation using singlet-triplet qubits was proposed by 
Levy [sj. Two-electron states in a double quantum 
dot (DQD) are characterized by the charge number in 
each dot labeled (m,n) and the total spin s = {S,T}. 
Then the basis set is composed of three singlet states 
\S±) = 0(2,0)5) ± |(0, 2)S)}/V2, \S) = |(1,1)5) and 
triplet states \T a ) = |(1, l)T a ) with a magnetic quantum 
number a — {0, ±1}. A single spin qubit is defined in one 
of the singlet-triplet subspaces, e.g., S-Tq subspace. The 
proposal has inspired different experimental [H, 0, [H, 0] 
and theoretical studies in recent years. 

In the S- To subspace, any of quantum superposition 
states can be written as 



cos-|S)+sm- e ^|T ) 



(1) 



and is mapped on the Bloch sphere in Fig. []Ja) . To build 
the quantum gates, it is necessary to manipulate and de- 
tect the Bloch angles 9 and (f>. However the conventional 
readout experiments using the Pauli spin blockade can- 
not detect the relative phase <fi, but can detect the total 
spin populations characterized by the angle 9 [H, 0, H, Q • 
The relative phase cf> is a fundamental element of the 
quantum mechanics in itself, and is essential in quantum 
algorithms such as Grover's database search problem 
Therefore we need to explore schemes for measuring the 
two Bloch angles in parallel. 

In this work, we propose a measurement scheme for 
the coherence phase between the entangled spin states, 
which utilizes adiabatic charge transfer. It is assumed 
that the DQD defined by metal surface gates is located 
on a two-dimensional electron gas (2DEG) in a semi- 
conductor. Then, we can control the electrical poten- 
tials in the two dots and the barrier potential separating 
them 0, [HI ■ We calculated the charge difference in 
the DQD when the quantum state in Eq. ([TJ) is tem- 
porally varied by the gate manipulation. It was shown 



that the expectation value of the transferred charge de- 
pends on the initial phase <f> due to spin-orbit couplings 
in the 2DEG. Thus, it becomes possible to identify the 
state vector on the Bloch sphere in the S-Tn subspac e by 
charge-sensing measurements at each dot 
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The system we consider is well described by the Hund- 
Mulliken model, in which doubly occupied states are 
taken into account. In this model, the orthonormal- 
ized single-particle state is defined as 3>l(R) = ( ( ( C L(R) — 
5 < /'r(l) )/ yl — 2Tg + g 2 where T is the overlap integral 
of the orbitals in the left and right dots </?l,r, and 
g = (1 — \/l — r 2 )/T [12]. We consider the spin-orbit 
interactions for the singlet-triplet coupling 13, 14, flij ]. 



Electron tunneling thus accompanies the spin preces- 
sion with respect to the vector P/2 = — ^^lI^J^r), 
in which the spin-orbit interaction is expressed by H so = 
Si = i.2^(k;) -Si with ki and Si being the wave vector and 
spin operator of z-th excess electron, respectively [l3| . It 
is convenient to take the spin quantization axis to be 
parallel to P/2. For this choice, the states \T±i) are de- 
coupled from the other states [13(, and will be omitted 
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FIG. 1: (color online) (a) Bloch sphere of the S-Tq spin sub- 
space. The angle 9 characterizes the total spin populations, 
whereas the angle <f) denotes the relative phase, (b) Scat- 
tering processes in the DQD; |5") and |Tb) couple with each 
other through the doubly occupied states. IS^c)) is a linear 
combination of |(2,0)S) and 1(0,2)5), which diagonalizes the 
doubly occupied states when we set T — P — virtually. 
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below. Thus, the Hamiltonian in {\S+), \S-), \S), \T )} 



basis can be written as 



H 



( U + X 

e 
T 

\ iP 



e T -iP\ 
U-X 





(2) 



in which U = (*l(r)|C|* L (r)> and X = (*l|C|*r> 
with * L(R) = $L(R)( r i)*L(R)('2) and C(ri,r 2 ) being a 
Coulomb interaction. T/2 and e are the single-particle 
interdot tunneling and the gate-controlled potential dif- 
ference between the dots, respectively. One can see that 
the singlet and the triplet (1,1) states couple with each 
other through the doubly occupied states [see Fig. [IJb)]. 
We assume that T and e are the accessible parameters 
by modulating the gate voltages. The matrix elements T 
and P are approximately proportional to T T3, Hj , and 
X to T 2 [HI in the tunneling regime. 

In Eq. |(3J), only the ground orbital of each dot is con- 
sidered. When one of the electrons is brought into the 
first excited orbital, the additional singlet and triplet 
states { | (0, 2)5') , | (2, 0)S'} , | (0, 2)T a ) , | (2, 0)T CT )} are pos- 
sible [8(. However, in the typical experiments, these 
states lie far ( > 0.4meV > T,P) above the \S±) 
state d, 0, @|. Thus we can disregard them as long 
as e < U. 

We performed a unitary transformation of \S±) so that, 
as |e| increases, one of them \S g ) energetically approaches 
the (1, 1) states while the other state \S C ) draws apart [see 
Fig- DIb)]. After adiabatic elimination of the higher state 
|<Se) [HI]) the effective Hamiltonian is given by 



Eg Ts'mr] —iPsmrj 

H e g — | Tsinri — 21 cos 2 rj i^- cos 2 r\ 
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iPsmr/ — i-g-cos rj 



(3) 
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where E, 



g(e) 



U T VX 2 + e 2 and tan?/ = (-X 



V 1 X 2 + s 2 )/e. Here we introduce the instantaneous 
eigenstates and energies of the time-dependent Hamil- 
tonian of Eq. such that H cS (t)\m(t)) = E m (t)\m{t)). 
Due to the presence of the spin-orbit couplings, the 
(1,1) charge state is expanded with the so-called bright 
state |B) = coseiS") + isine|r ) and dark state |D) = 
sin BIS") - icos6|r ) with tan 6 = P/T 18]. Therefore 
one sees that the instantaneous eigenstates are |D) and 



|+) = cos$|Sg) +sin$|B) 
|-) = -sin^l^g) +cos$|B) 



in which 



tan $ = 



E„ — E„ 



E 2 + 4(T 2 + P 2 ) sin 2 ri 



v/4(T 2 + P 2 )sinr? 
(T 2 + P 2 ) 



■ cos 2 rj. 



(4) 

(5) 
(6) 



The dark state |D) is free from the double occupancy 
state, and is not affected by the sweep of the bias poten- 
tial. In addition, since the spin-orbit coupling is always 
weak compared with the hopping energy, the variation 
of the mixing angle 8 by the center gate is considerably 
small. On the other hand the angle $ changes from 7r/2 
to with increasing e, which indicates that there is an 
avoided crossing of the |B) and \S g ) states. The corre- 
sponding eigenenergies are respectively 



E D =0, 



E± = 



E g =p \ E 2 + 4(T 2 + P 2 ) sin 2 r? 



(T 2 + P 2 ) 



(7) 



cos 2 ry. 



We calculate the charge difference in the DQD after 
varying the system parameters adiabatically. The ex- 
pected value of the charge difference is defined by 



Q(r) = / dt(rP(t)\r\iP(t)), 
Jo 



(8) 



in which V is 3 x 3 the "current" operator. The cur- 
rent operator is obtained by mapping I = (i/h)[H, p z ] on 
{|D),|+),|-)} basis, where p z = 2e[|(2,0)5)<(2,0)5| - 
|(0,2)S')((0,2)5 , |] and [ , ] denotes commutation. Be- 
cause the dark state |D) is decoupled from the double 
occupancy state, it does not contribute to the current. 
Then the detection of Q(t) corresponds to the projec- 
tion measurement in the bright state |B). Note that in 
the Pauli spin blockade measurement the projection axis 
is \S) [H, 0> M [H • Since the bright state lies on the sphere 
with a well-defined azimuthal angle of </>b = 7r/2, the 
charge difference can capture the initial relative phase. 

With use of the instantaneous eigenstates, one can ex- 
pand a two-electron state as 



\m)= E « 

m={D,±} 



(t)e*»|m;q), 



(9) 



where Cm(t) = — L dt'E m {t')/h is the usual dynamical 
phase. The coefficient c m is varied with respect to the 
set of gate-controlled paramete rs qft) = {T{t),e(t)}, and 
obeys the differential equation 



rli ~ ~ 2^, 



{m;q\ — \n;q) 



(10) 



The time variation of 5c m (t) = c m (0) — c m (t) represents 
the nonadiabatic level transition. However, as far as the 
adiabatic condition |(m; q\h(d/dt)\n; q) /(E m — E n )\ -C 1 
is satisfied, it is negligible and may be dropped. 

The adiabatically pumped charge difference is obtained 
by substituting Eqs. (|9ll0p into Eq. ©, and taking the 
zero-order terms in terms of Sc m (t). Here we assume 
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FIG. 2: (color online) Sequence of gate voltage manipulations 
under consideration. The upper panel shows the potential 
bias between the left and the right dots. The overlap of the 
orbitals in the two dots V is moved up and down by the gate- 
controlled potential barrier. 



that the manipulation time r is much longer than the 
period of the unitary time evolution e l< > m . As a result, 
we find that the charge difference consists of two parts; 
Q(t) = Q° + Q\ One is 



Q° = 2Hm V |c m (0)| 2 f dt 

„ i_ Jo 



T ^d,, (m\P\n) 



[n\ — \m) 



dt E. m — E n 



(ii) 



The other is the interference part, which includes rapid 
inter-level oscillation in the integrand; 

Q 1 = 2fdm V c m (0)<(0) / dte^-^K mn , (12) 
™ „j™ Jo 



where 



K mn (t) = —{n\I'\m) 



E 



, . d , , (l\F\m) . „. 



When an eigenstate is prepared as an initial state, i.e., 
Cm(0) = 1 for a certain m and Cm'y^m (0)_= 0, Q(t) re- 
duces to the result in the previous work 1£J, |2fJ • 

Hereafter we calculate the charge difference Q(t) for 
specific manipulation sequence. The gate voltages are ini- 
tially adjusted so that no bias potential is present, and 
the barrier potential is so high that the system stays al- 
most in the prepared state. The gate control under con- 
sideration is presented in Fig. O The sequence consists 
of three parts; we (i) lower the interdot potential barrier 
in order to increase the overlap integral T, (ii) tilt the 
electric potential until e = ei > to be E s (s{,Tq) = 0, 
and (iii) raise the barrier potential height again. 

Before proceeding to the calculation, we review the adi- 
abatic conditions. First, no charge transfer occurs dur- 
ing the process (i) because the doubly occupied state is 
entirely decoupled. In the second process (ii), the gate 
sweep affects only the mixing angle $. The dark state |D) 
is thus left unperturbed. Then, the nonadiabatic contri- 
bution Sc m is roughly proportional to the Landau- Zencr 



(LZ) transition rate 2l| between the doubly occupied 
state \S g ) and the bright state |B). At the avoided cross- 
ing, the LZ transition rate is estimated as 

en) „„r <n + Pi)^ 



V 



expi 



h\e\ 



(14) 



where i = de/dt, T 0(1 ) = T(r 0(1) ) and P 0(1) = P(F 0(1) ). 
In the last process (iii), lowering T changes both the 
mixing angles and $, but generally dQ/dt <C d^/dt. 
Thus, assuming that the applied bias is so large that 
rjf(t) = rj(ef, T) ~ 7r/4, the LZ transition rate is 



P 



(iii) 



exp 



j- TrE c (e { ) 2 \ 
V" ft|A| J 



(15) 



where A = dy/T 2 + P 2 /dt is approximately proportional 
to dY/dt. It should be noted that the spin-orbit inter- 
action appears in parallel with the interdot tunneling as 
long as the dark state is robust against the gate con- 
trol. The coherent oscillation between the singlet and the 
triplet states holds when the adiabatic condition for the 
interdot tunneling T/2 is satisfied. This is a definite dif- 
ference between the previous work using a field gradient 
for singlet-triplet mixing [1,0,0,0,0] and ours. Besides, 
the single level spacing for the individual electrons is of 
the order of lmeV \JT 2 + P 2 . Therefore, the nona- 
diabatic transition to the excited orbitals is negligible as 
long as , p( m ) <C 1 is satisfied. 

The target state we focus on is arbitrary superposi- 
tion state of \S) and |To) as shown in Eq. ([T]). At the 
initial condition t = 0, it does not contain the excited 
state |— ), i.e., c_(0) = 0. Therefore one can see that the 
interference part becomes 

, [ T e iC +sin$ d , s 

Q 1 cx / dt ; =~tQ- (16) 
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E 2 + 4(T 2 + P 2 ) sm 2 ry 



As is mentioned above, the integrand is negligibly small 
and rapidly oscillating. Thus all we have to calculate 
is Q°, which is proportional to |c + (0)| 2 . Within the first 
order of P/T, we obtain the adiabatically pumped charge 
difference as 



QW = -|{(i 



p 

cos 0) + 2— sin( 
To 



b). (17) 



The first term in the right-hand side describes the Pauli 
spin blockade d, 0, 0, H3| ■ It should be noted that the 
charge difference oscillates with respect to the relative 
phase 0. The oscillating term reflects the imaginary part 
of \ip), and originates from the projection axis |B) which 
is slightly tilted from 15). In addition, the oscillation 
does not appear when the initial state is a pure singlet 
(9 = 0) or triplet state (9 = n). 

For a clean GaAs/AlGaAs 2DEG confined in 10 — 
lOOnm long, the spin-orbit interaction energy is esti- 
mated to be (H S o) = 1CP 2 ~ 10 _1 meV from magne- 
toresistance data [23j. On the other hand the confine- 
ment energy is ~ lmeV for a quantum dot with a 30nm 
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side [12|]. In that case, it is possible to experimentally 
achieve the condition in which the oscillation amplitude 
is - 10% of Q(t,</> = 0) [H. Therefore the repetitive 
experiments can reveal the relative phase as well as the 
total spin population of the prepared state. 

So far, we have neglected the effect of nuclear spins 
in the semiconductor. The hyperfine fields due to the 
nuclei (Overhauser field) h.L(R) couple with the electron 
spins as 14f = h • (Si + S 2 ) + 5h ■ (Si — S2), where 
h = (h L + h R )/2 and 5h = (h L - h R )/2. Then the 
average of the Overhauser fields h rotates the subspace 
of the three spin triplet states, while the inhomogeneity 
Sh mixes \S) with \T a )s ^M,^- 

The Overhauser fields can disturb the electron spin 
state during the adiabatic gate control. However the hy- 
perfine coupling between \S) and \T a ) does not undergo 
virtual double occupancy, and the adiabatic condition for 
5h is different from that for T/2 and P/2. Thus, we can 
separate off the effect of the Overhauser field using the 
technique called "rapid adiabatic passage" , in which the 
sweep of the bias is adiabatic for the electron tunneling 
but is nonadiabatic for the hyperfine couplings 0, (2?J ■ In 
a quantum dot containing unpolarized N = 10 5 nuclear 
spins, the root mean square of the Overhauser field is 
|(h L ( R )) rms | ~ 10 _4 meV. The required length of the ma- 
nipulation sequence r is a few /is or shorter for interdot 
tunneling coupling Ti/2 ~ 10 _2 meV. This condition has 
been achieved in a couple of experiments 0, • 

In summary, we propose an adiabatic charge transfer 
in a gate-defined DQD as an indicator of singlet-triplet 
quantum superposition on a S-Tq Bloch sphere. After 
the gate manipulations, the transferred charge number 
is found to oscillate with respect to two Bloch angles in 
the initially prepared superposition state. The oscilla- 
tion can be observed in an ensemble average of charge- 
sensing measurements in each dot with quantum point 
contacts 0, 0, EH • Recently Kosaka et al. demonstrated 
the quantum coherence transfer from light polarization 
to electron spin polarization in a quantum well 28|, 2^ |. 
By applying this method, it becomes possible to prepare 
arbitrary S-Tq superposition states in DQD. The present 
scheme can help to check whether a system is indeed pre- 
pared in the desired state. 
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